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oo ! Abstract 

Field theoretical renormalization group methods are applied to the Obukhov- 

^ Kraichnan model of a passive scalar advected by the Gaussian velocity field 

& . 

^ c~| with the covariance 

O 

(v(t,x)v(t',x)) - (v(t,x)v(t',x')) oc 5(t-t')\x-x'\ £ . 

Inertial range anomalous scaling for the structure functions and various pair 
correlators is established as a consequence of the existence in the correspond- 
ing operator product expansions of "dangerous" composite operators [powers 
of the local dissipation rate], whose negative critical dimensions determine 
anomalous exponents. The main technical result is the calculation of the 
anomalous exponents in the order e 2 of the e expansion. Generalization of 

the results obtained to the case of a "slow" velocity field is also presented. 
PACS number(s): 47.10.+g, 47.25.Cg, 05.40.+j 
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I. INTRODUCTION 



One of the main problems in the modern theory of fully developed turbulence is to 
verify the basic principles of the Kolmogorov-Obukhov (KO) phenomenological theory 
within the framework of a microscopic model and to investigate deviations from this theory, 
provided they exist. 

In particular, one is interested in the single-time "structure functions" 

S„(r) = <[0(x)-0(x')P>, r=|x-x'| (1.1) 

in the inertial range. Here 9{x) = 9(t, x) can be the component of the velocity field directed 
along the vector x — x', or the scalar field in the problem of turbulent advection; the brackets 
(. . .) denote the ensemble averaging, and the time argument common to all the quantities is 
omitted in ( |1.1|) and analogous formulas below. The inertial range (or the convective range 
in the problem of turbulent advection) is defined by the inequalities I << r << L, where 
I = A -1 is the internal (viscous) scale and L = M _1 is the external (integral) scale. 

According to the KO theory, the functions ( |1.1|) in the inertial range are independent 
of both the external and internal scales (the First and the Second Kolmogorov hypotheses, 
respectively) and are determined by the only parameter e, the mean dissipation rate, see 
e.g. p]|3]. Dimensionality considerations then determine the functions ( |1 . 1| ) , apart from 
numerical coefficients, in the form 

S n (r) = const (er) n/3 . (1.2) 

Both experimental and theoretical evidence is known in favor of some deviation from the 
predictions of the KO theory, see |],[3|-[{j] . For the structure functions Ql.l]) , these deviations 
are phenomenologically written in the form (in contradiction with the First Kolmogorov 
hypothesis) 

S n (r) = const (tr) n ' z f n (Mr), (1.3) 

where the "scaling functions" f n (Mr) are supposed to have power-like behaviour in the 
asymptotic region Mr << 1, 

fn(Mr) ~ const (Mr) Qn . (1.4) 

The singular dependence of the structure functions on M for M — > and nonlinearity 
of the exponents q n in n are usually referred to as "anomalous scaling", and in theoretical 
models they are explained by strongly developed fluctuations of the local dissipation rate 
( "intermittency" ) , see discussion in IJ, Chap. 25. Within the framework of numerous models, 
the anomalous exponents q n are related to the statistics of the local dissipation rate or to the 
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dimensionality of fractal structures formed by small-scale vortices in the dissipative range, see 
for a review 0]. As a rule, these models are only loosely related to underlying microscopic 
models, and therefore some doubt remains about the universality of representations like 
3D , ( |1.4J) and the very existence of deviations from the KO theory. 

An effective method of studying self-similar scaling behaviour is that of the renormaliza- 
tion group (RG), see ||||. It was successfully applied in the theory of critical phenomena to 
explain the origin of critical scaling and to calculate universal quantities (critical dimensions 
and scaling functions) in the form of series in the formal small parameter e = 4 — d, where 



d is the space dimensionality, see ||10|| . This technique is also fully applicable to the theory 



of turbulence ||11|| , see also the review paper |L2| and references therein. 



In the statistical theory of turbulence, the microscopic model is usually taken to be the 
stochastic Navier-Stokes equation (SNS) with an external random force which imitates the 
injection of energy by large-scale modes, see e.g. [|iTJ^|r2f . The role of the RG expansion 
parameter is played by the exponent e entering into the random force correlator. For the 
structure functions (|1 . 1| ) , the RG method allows proof of the Second Kolmogorov hypoth- 
esis (independence of the viscous length) for a wide variety of realistic random forces, see 
fT^i4|. This is equivalent to the representation ( |1.3| ), in which the form of the scaling func- 



tion f n (Mr) is arbitrary. We note that the representations (|1.3j ) for any functions f n (Mr) 
imply the existence of scaling in the infrared (IR) region L,r >> I with definite "critical 
dimensions" 

A M = -A r = 1, A[3 n ] = -n/3 (1.5) 



for any fixed value of e. This means that the structure functions ( |1.3| ) scale as S n — > X A ^- Sn ^S n 
upon the substitution M — > X Am M, r — > X Ar r [in general, the exponent A[S n ] is replaced by 
the critical dimension of the corresponding correlation function; this dimension is calculated 
within the e expansion]. It is this IR scaling that is analogous to the critical scaling in the 
theory of critical phenomena, where the role of an external scale is played by the correlation 
length L = r c , see ||10|| . In the RG approach, the critical dimensions like (|1.5|) arise as 
coefficients in the RG equations. The form of the scaling functions is not determined by the 



RG equations themselves, and therefore the anomalous exponents in (1A) are not related to 
the critical dimensions of the functions ( |1.1| ). 

As in the case of critical phenomena [|IIJ] , the dependence of the scaling functions on the 
argument Mr in the region Mr << 1 is studied here using the well known operator product 
expansion (OPE, or SDE=short distance expansion), see [12-17]. According to the OPE, 
the scaling functions in expressions like ( |1.3|) are represented in the form 



/(Mr) = J2 C F(Mr) AF , Mr -> 0. (1.6) 
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Here Cf are coefficients regular in M, the summation is implied over all possible composite 
operators F allowed by the symmetry of the left-hand side, and Ap are their critical di- 
mensions. In particular, only scalar Galilean invariant operators with nonzero mean value 
contribute to the OPE for the structure functions ( |1.1|) , see flT^JTSU . 

In the theory of critical phenomena, all the nontrivial composite operators have positive 
critical dimensions, > 0, and the leading term in ( |1.6| ) is determined by the simplest 
operator F = 1 with Ap = 0, i.e., the function /(Mr) is finite as M = r^ 1 — > 0, see 
|T0| . It has long been realized (see also papers [12,14-17]) that the singular behaviour 
of the scaling functions in the SNS model for M — > is related to the existence in the 
model of composite operators with negative critical dimensions. These operators, called 



"dangerous" in |C| , should not be confused with IR relevant operators considered in fT8| . |19 
whose existence can lead to the violation of the scaling regime. 

Dangerous composite operators in the SNS model occur only for finite values of the RG 
expansion parameter e, and within the e expansion it is impossible to decide whether or not 
a given operator is dangerous, provided its critical dimension is not found exactly using the 
Schwinger-type functional equations or the Galilean symmetry, see [12,15-17,19-21]. More- 
over, dangerous operators enter into the operator product expansions in the form of infinite 
families with the spectrum of critical dimensions unbounded from below, and the analysis 
of the small M behaviour implies the summation of their contributions. Such a summation 
for the case of different-time correlators, first accomplished in |l3j (see also |l2|p!5| ), estab- 
lishes the substantial M dependence of the correlators and their superexponential decay as 
the time differences increase. Of course, this effect is well-known and has a simple physical 
interpretation as the transport of turbulent eddies as a whole by the large-scale modes, see 
|22| . Therefore, this effect is not ignored within the correct RG formalism, as it is sometimes 



believed |23|J24j] , provided the RG is used beyond the e expansion and is combined with the 
OPE technique. 

The analysis of the M dependence of Galilean invariant quantities, such as the single- 
time structure functions Ql.lp , requires the explicit construction of all dangerous invariant 
scalar operators, exact calculation of their critical dimensions, and summation of their con- 
tributions in the corresponding operator product expansions ( |1.6[ ). This is clearly not a 
simple problem and it requires considerable improvement of the present technique. 

In view of the difficulties encountered by the RG approach to the SNS model it is tempt- 
ing to apply the formalism to simpler models, which exhibit some of the features of real 
turbulent flows, but are easier to study. An interesting example is provided by the well- 
known Heisenberg model, see Chap. 17, whose exact solution has recently been rederived 
using the RG method [25|. Unfortunately, the Heisenberg model does not involve higher- 
order functions and is therefore not suitable for studying anomalous scaling. 

Recently, much attention has been attracted by a simple model of the passive advection 
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of a scalar quantity by a Gaussian velocity field, introduced by Obukhov p6 and Kraichnan 



27], see the papers [28-40] and references therein. It turns out, that the structure functions 
of the scalar field in this model exhibit anomalous scaling behaviour analogous to ( |1.3| ), fll.4| ), 
and the corresponding anomalous exponents can be calculated explicitly within expansions 
in certain small parameters, [31-34]. 

The advection of a passive scalar field 8(x) = 6(t,x) is described by the stochastic 
equation 

V t e = u A6 + f, Vt^dt + Vtd,, (1.7) 

where dt = d/dt, di = d/dxi, vo is the molecular diffusivity coefficient, A is the Laplace 
operator, v(x) is the transverse (owing to the incompressibility) velocity field, and / = f(x) 
is a Gaussian scalar noise with zero mean and correlator 

(f(x)f(x')) = 6{t-t')C{Mr), r=|x-x'|. (1.8) 

The parameter L = M~ l is an integral scale related to the scalar noise, and C(Mr) is some 
function finite as L —>■ oo. Without loss of generality, we take C(0) = 1 (the dimensional 
coefficient in ( |1.8| ) can be absorbed by appropriate rescaling of the field 9 and noise /). 



In a more realistic formulation, the field v(x) satisfies the stochastic Navier-Stokes equa- 
tion, see e.g. [4~I|. Following [26-34], we shall consider a simplified model in which v(x) obeys 



a Gaussian distribution with zero average and correlator 

(v^vjix')) = D 6 -^-p- J dki^-(k) (k 2 + m 2 )- d ^ 2 exp[ik ■ (x - x% (1.9) 

where -£\j(k) = 5ij — kikj/k 2 is the transverse projector, k = |k|, Dq is an amplitude factor, 
1/m is another integral scale, and d is the dimensionality of the x space; < e < 2 is a 



parameter with the real ( "Kolmogorov" ) value e = 4/3; e = 2 — 7 in the notation of |31"| , |3~2" 
e = k in the notation of |33|], and e = \ m the notation of [34fl. The relations 



g = D /u = A e (1.10) 

define the coupling constant go (expansion parameter in the perturbation theory, see Sec.|T|) 
and the characteristic UV momentum scale A (A ~ l/r^ in the notation of pi|j5^ ]). 

In the model ( |1.7| )-( |i~9"D , the odd multipoint correlation functions of the scalar field 
vanish, while the even single-time functions satisfy linear partial differential equations. The 
solution for the pair correlator is obtained explicitly; it shows that the structure function 



S2 is finite for M,m = | 27fl . The higher-order correlators are not found explicitly, but 
their asymptotic behaviour for M — > can be extracted from the analysis of the nontrivial 
zero modes of the corresponding differential operators in the limit 1/d — > 0,0] or e — > 
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33|,|34|. It was shown that the structure functions are finite for m = 0, and in the convective 



range A >> 1/r >> M they have the form (up to the notation): 

S 2n {r) = <[0(x) - 0(x')] 2n ) oc D n r n ^ (Mr ) An , (1.11) 



with negative anomalous exponents A„, whose first terms of the expansion in 1/d []3T|J32 
and e fl33lj34|[ have the form 

A n = -2n(n - l)e/(d + 2) + 0(^ 2 ) = -2n(n - l)e/d + 0(l/d 2 ). (1.12) 



We note that A n = — p 2n in the notation of [33], A 2 = —A in the notation of [31 1, and 
A n = -A 2n in H. 

Another quantity of interest is the local dissipation rate, E(x) = h> {di6(x)di6(x)}. The 



single-time correlation functions of its powers in the convective range have the form p2| , |34 

(E n (x) E m (x')) oc (Ar)- A "- Am (Mr) A " +m (1.13) 

where A is defined in ( |1.10| ) and the exponents A n are the same as in Eq. ( |1 . 1 1| ) . 

In this paper, we apply the RG and OPE approach developed in [12-17,19-21,41] to 
the model ( |1.7] )-( rO| ). We show that the RG explains the origin of anomalous scaling and 
allows the anomalous exponents to be calculated in the form of a series in the parameter 
e entering into the correlator fll.Sf ). Therefore, this parameter plays in the RG approach 
the role analogous to that played by the parameter e = 4 — d in the RG theory of critical 
phenomena, and the results obtained in |33],[34| can be interpreted as the first terms of the 
corresponding e expansions. 

We have calculated the exponents A n in the second order of the e expansion for an 
arbitrary value of the space dimensionality d (they are given in Sec.|lll|). In particular, we 
have obtained 

n(n — l)e n(n — l)e 2 r , , . , .. , .,, q . 

A n = 1 - ; + 1 - J 2n 19 - 66 In 4/3 +5-24 In 4/3 + O (e 3 ) = 

I 16 

= -n(n - l)e/2 + e 2 n(n - l)[0.00162n - 0.11902] + 0(e 3 ) (1.14) 



for d = 2 and 

2n(n-l)e 2n(n - l)e 2 

A n = — H i — x 

5 875 

x 2n(2557rV3 - 1384) + 345^^ - 1884] + 0(e 3 ) = 
= -2n(n - l)e/5 + e 2 n(n - l)[0.01626n - 0.01535] + 0(e 3 ) (1.15) 

for d = 3. 

Representations like ( |1 . 1 1| ) , ( |1.13| ) can also be derived using the RG and OPE for any 
correlation function; as an additional example, we consider the single-time pair correlators 
involving second-rank irreducible tensors of the form 
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a i ed j e[d.ed.e) n - 2 - 5 %J [d s ed s d} n /d, 



the special case n = 2 was considered previously in p2| . 

From the RG viewpoints, the model ( |1.7| )-( |L9| ) is simpler than the "true" SNS model 
in at least two respects. First, dangerous operators in the model ( |1.7|) -( |L9D exist even for 
asymptotically small values of e and therefore they can be identified within the e expansion. 
These are the local dissipation rate E(x) and all its powers, their critical dimensions being 
A n . Second, only finite number of these operators contribute to the operator product expan- 



sion (1^6) for any given correlation function, and the additional resummation of the series 
discussed above, is not required here: the leading term of the asymptotic behaviour of 
the function f(Mr) for Mr — > is simply given by the contribution of the "most dangerous" 
operator, i.e., that having the smallest value of A n . 



It was suggested in [33[] that the anomalous exponents describing the asymptotic be- 
haviour for Mr << 1 can be related to certain IR relevant interactions in a hypothetical RG 
approach, see also discussion in [32]. This is not exactly so: in fact, these exponents in the 



model (|1.7| )- (|1.9| ) are related to the critical dimensions of certain "dangerous" composite 
operators in the operator product expansions for corresponding Green functions. It should 
be stressed, that in the SNS model the behaviour at Mr — > will probably be determined 
by additional resummations of the operator product expansions, which can lead to more 
exotic behaviour than the simple power-like one in ( |1 . 1 1|) , ( p.. 13 ). 



The plan of our paper is the following. In Sec.||, we discuss the field theoretical formu- 
lation and UV renormalization of the model (|1.7| )- (|OD and derive the corresponding RG 
equations with exactly known RG functions (the (3 function and the anomalous dimension). 
These equations have an IR stable fixed point, which establishes the existence of IR scaling 
with exactly known critical dimensions of the basic fields and parameters of the model. In 
Sec.p|, we discuss the renormalization of various composite operators in the model (|1.7|) - 
(|1.9|). In particular, we present the second-order result for the e expansion of the quantity 



A n , the critical dimension associated with the operator E n (x). In Sec.|IV|, it is explained 
using the OPE that, at the same time, A n play the part of the anomalous exponents in 
relations like (|1.11| ), ( |1.13j ). The results obtained are discussed in Sec.|V|. We also briefly 



mention there one of the possible modifications of the model ( |1.7| )-( |i"D| ), that of a "slow" 
velocity field, in which the correlator (|1.9| ) contains no 5 function in time. The RG analysis 
shows that the Green functions in this model exhibit an anomalous scaling behaviour, and 
the corresponding anomalous exponents are calculated in the form of series in e + 2. 
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II. FIELD THEORETICAL FORMULATION, UV RENORMALIZATION, AND 

RG EQUATIONS 



The stochastic problem ( |1.7p -( |OD is equivalent to the field theoretical model of the set 
of three fields $ = {6, 6', v} with action functional 

s-($) = e'D e e'/2 + & \-d t e - (yd)e + v A6] - v d; 1 v/2. (2.1) 

The first four terms in ( [2.1| ) represent the Martin-Siggia-Rose-type action [42-45] for the 
stochastic problem ( |1.7|) , (|1.8| ) at fixed v, and the last term represents the Gaussian averaging 
over v. Here Do and D v are the correlators ( |i.8p and respectively, the required 



integrations over x = (t, x) and summations over the vector indices are understood. 

The formulation (|2.1|) means that statistical averages of random quantities in stochastic 
problem ( |1.7| )-( |1~9|) coincide with functional averages with the weight expS'($), therefore 
generating functionals of total (G(A)) and connected (W (A)) Green functions of the problem 
are represented by the functional integral 

G(A) = exp W(A) = J £>$ exp[5($) + A$] (2.2) 

with arbitrary sources A = A 9 , A 9 ' , A v in the linear form 

A$ = J dx[A e {x)6{x) + A e \x)6\x) + AJ{x)vi{x)}. 



The model ( ]2.1|) corresponds to a standard Feynman diagrammatic technique with the 
triple vertex —6'(yd)6 and bare propagators (in the momentum-frequency representation) 

(ee')o = (e'6)* = {-iu + u k 2 )-\ 

)o = C(k) (u 2 + vlk'Y\ 

(O'e')o = 0, (2.3) 



where C(k) is the Fourier transform of the function C(Mr) in ( |1.8| ) and the bare propagator 
(vv)o is given by Eq. (|1.9|). The role of the coupling constant in the perturbation theory is 
played by the parameter go defined in ( |1.10|) . 

It is well known that the analysis of UV divergences is based on the analysis of canonical 
dimensions. Dynamical models of the type (|2.1|) , in contrast to static models, are two-scale, 
i.e., to each quantity F (a field or a parameter in the action functional) one can assign 
two independent canonical dimensions, the momentum dimension d k F and the frequency 
dimension d%, determined from the natural normalization conditions 

d\ = -di = = < = 0,d* = d\ = 0, dl = -d% = 1, 



S 



and from the requirement that each term of the action functional be dimensionless (with 



respect to the momentum and frequency dimensions separately), see |jO|,20|. Then, based 
on dp and d F , one can introduce the "summed" (total) canonical dimension dp = d k F + 2d F 
(in the free theory, dt oc A). 

The dimensions for the model ( |2.1| ) are given in Table |, including renormalized param- 
eters, which will be considered later on. 

From Table | it follows that the model is logarithmic (the coupling constant go is dimen- 
sionless) at e = 0, and the UV divergences have the form of the poles in e in the Green 
functions. 

The total dimension dp plays in the theory of renormalization of dynamical models the 
same role as does the conventional (momentum) dimension in static problems. The canonical 
dimensions of an arbitrary 1-irreducible Green function T = ($ . . .$)i_; r are given by the 
relations 



d^= d - N^d®, 
d% = l- N$d%, 

d r = 4 + 2d^ = d + 2 - JV$d$, (2.4) 

where = {Ng, Ngi, N v } are the numbers of corresponding fields entering into the function 
T, and the summation over all types of the fields is implied. The total dimension dr is the 
formal index of the UV divergence. Superficial UV divergences, whose removal requires 
counterterms, can be present only in those functions T for which dr is a nonnegative integer, 
see e.g. 

Analysis of divergences should be based on the following auxiliary considerations: 

(1) From the explicit form of the vertex and bare propagators in the model ( |2.1|) it follows 
that Nqi — Ng = 2N for any 1-irreducible Green function, where N > is the total number 
of the bare propagators (99)o entering into the function (obviously, no diagrams with N < 
can be constructed). Therefore, the difference Ngi — Ng is an even nonnegative integer for 
any nonvanishing function. 

(2) If for some reason a number of external momenta occur as an overall factor in all the 
diagrams of a given Green function, the real index of divergence d' T is smaller than dr by 
the corresponding number of unities (the Green function requires counterterms only if d' r is 
a nonnegative integer). 

In the model (|2.1| ), the derivative d at the vertex 9'{wd)9 can be moved onto the field 
9' by virtue of the transversality of the field v. Therefore, in any 1-irreducible diagram 
it is always possible to move the derivative onto any of the external "tails" 9 or 9' , which 
decreases the real index of divergence: d' r = dr — Ng — Ngi. The fields 9, 9' enter into the 
counterterms only in the form of derivatives d9, 89'. 
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From the dimensions in Table | we find dr = d + 2 — N v + Ng — (d + l)N$i and 
d' r = (d + 2)(1 — iVfl/) — N v . From these expressions it follows that for any d, superfi- 
cial divergences can exist only in the 1-irreducible functions (O'O . . .9) with Ng/ = 1 and 
arbitrary value of Ng, for which dr = 2, d' r = 0. However, all the functions with Ng > Ng/ 
vanish (see above) and obviously do not require counterterms. We are left with the only 
superficially divergent function (O'O); the corresponding counterterm must contain two sym- 
bols d and is therefore reduced to 8' AO. Inclusion of this counterterm is reproduced by the 



multiplicative renormalization of the parameters go , Uq in the action functional (2.1) with 
the only independent renormalization constant Z v : 

v = vZ u , g = gn 6 Z g , Z g = Z~ x . (2.5) 

Here // is the renormalization mass in the minimal subtraction scheme (MS), which we always 
use in what follows, g and v are renormalized analogues of the bare parameters go and vq, 
and Z = Z(g, e, d) are the renormalization constants. Their relation in fl2.5|) results from the 
absence of renormalization of the contribution with Do in Q2.1p , so that -Do = ^o^o — 9^ £v i 
see ( P-.lOp . No renormalization of the fields and "masses" is required, i.e., Z§ — 1 for all $ 



and m = m, M = M, Z m = Zm = 1- 

Since the fields are not renormalized, their renormalized Green functions W R coincide 
with the corresponding unrenormalized functions W = ($...$) (for definiteness, we discuss 
the connected functions); the only difference is in the choice of variables and in the form of 
perturbation theory (in g instead of g ): 

W R {g ) v^ ) ...) = W{go,voT..) (2.6) 

(the dots stand for other arguments like coordinates and momenta). We use to denote 
the differential operator [id^ for fixed bare parameters go, uq and operate on both sides of 
Eq. ( p.6| ) with it. This gives the basic differential RG equation: 

V RG W R (g,u,fi,...) = 0, 
V RG = V^ + (3(g)d g - lv {g)V v , (2.7) 

where we have written T> s = sd s for any variable s, and the RG functions (the (3 function 
and the anomalous dimension 7) are defined as 

7^=^111^, (3{g) = V ll g = g[-e + lu ]. (2.8) 

The relation between (3 and 7 results from the definitions and the last relation in (|2.5|) . In 
general, if some quantity G is renormalized multiplicatively, G = Z G G R , it satisfies the RG 
equation of the form 
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[V RG + j G (g)} G R = 0, lG (g) = In Z G 



(2.9) 



with the operator T>rg from 

Explicit calculation of the constant Z v in the model (|2.1|) in the one-loop approximation 
gives: 



1 - 



g(d-l) C d 
2de 



(2.10) 



where we have written C d = S d / {2ix) d and S d = 2n d ^ 2 /T(d/2) is the surface area of the unit 
sphere in rf-dimensional space. 

The one-loop approximation Q2.10 ) for the constant Z v is in fact exact, i.e., it has no 
corrections of order g 2 , g 3 , and so on. Indeed, from the explicit form of the vertex and 
the bare propagators (|1.9| ), ( |2.3| ) it follows that any multiloop diagram of the 1-irreducible 
function (6'6) contains effectively a closed circuit of retarded propagators (99')q and therefore 
vanishes (it is also important here that the propagator (vv)o in ( |1.9j ) is proportional to the 
S function in time). 

From the definitions ( |2.8|) using Eq.( |2.1U[ ) we find exact expressions for the basic RG 
functions: 



g(d- l)C d M v 

lAg) = ' P( g ' = 9 



—e 



g{d-l)C d 
2d 



(2-11) 



From ( |2.11| ) it follows that an IR-attractive fixed point 

2de 



9* 



C d (d-1] 



(2-12) 



of the RG equations \j3(g*) = 0, /3'(g*) = e > 0] exists in the physical region g > for all 
e > 0. The value of 7^(5') at the fixed point is also found exactly: 



it = lu{g*) = e, 



(2-13) 



without corrections of order e 2 , e 3 , and so on. 

It is well known that the leading term of the IR asymptotic expression of any renormalized 
quantity G R , for which the RG equation of the form (|2.9|) is valid, satisfies the same equation 
with the substitution g —> g*, where g* is the IR stable fixed point: 



Canonical scale invariance is expressed by the relations: 



d% 



G 



0. 



G 



0. 



(2-14) 



(2-15) 
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in which a = {t, x, /i, is, m, M, g} is the set of all arguments of G R (t, x is the set of all times 
and coordinates), and d k and d w are the canonical dimensions of G R and a. Substituting 



the needed dimensions from Table | into Q2.15), we obtain: 



V u + V, 



V 



2V V -V^-d k 



G 



G R = 0, 



d G ) G R = 0. 



(2.16) 



Each of the equations (|2.14 )— ( 2.16 ) describes scaling with dilatation of all variables, the 
derivatives with respect to which enter into the differential operator. We are interested in 
scaling with dilatation of t, x, and "masses" M, m for fixed /i, is, and g, and it is necessary to 
exclude the corresponding derivatives T> a by a combination of the available equations. After 
eliminating and T> v from ( |2.14| ) and ( |2.16| ) we obtain the desired equation of critical IR 
scaling for the model (|2.1|): 



[-X? x + A t V t + A m V m + A M V M - A G ] G 



R 







with the coefficients 



At 



-A, 



-2 + e, A T 



My 



and 



A[G] = A G = d k G + A u d% + >y* G , 



(2.17) 



(2.18) 



(2.19) 



which are the corresponding critical dimensions. In particular, for any correlation function 
G R = W R = ($...$) of the fields $ we have A G = N^A^,, with the summation over all 
fields $ entering into G R , and for A$ using the data from Table [I] and the exact value of 
7* = £ we obtain from ( |2.19| ) the following exact expressions: 



A v 



1-e, A e = -l+e/2, A e , = d+ l-e/2 



(2.20) 



for 



(we recall that the fields in the model fl2.1p are not renormalized and therefore 7$ 
all $). 

To avoid misunderstanding, we again emphasize the fact that the RG equation (|2.17|) 
describes IR scaling, i.e., the statement which is equivalent to the critical scaling in the 
theory of critical phenomena, see also discussion in [|l2| , |i"3f . In this scaling the variables M, 
m are also IR relevant, i.e., they are dilated in scale transformations. In other words, the 
expression (|2.17 ) describes the asymptotic behaviour as Ar >> 1 for any fixed values of 
Mr and mr, with the UV scale A defined in (1.10). The solution of the set of equations 
(|2.16|) , ( 2.17|) can be found only up to some unknown function of all the first integrals, 



including those of the form mr, Mr, where r 



x,-| is some coordinate difference. 
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The anomalous exponents in the expressions ( 1.11J) , Ql-13|) describe the behaviour of the 
corresponding correlation functions for Mr — > 0, and therefore are not related to the their 
critical dimensions ( [2.19| ). 

In terms of unrenormalized variables, the solution of the set of equations (|2.15|) , ([2.171) 
for the example of a single-time pair correlation function can be found, apart from numerical 



factor, in the form (for more details, see e.g. fll2|): 



GocG 1 



f(Mr, mr), 



(2.21) 



with certain, as yet unknown, scaling function /. 

In what follows we limit ourselves to the correlation functions of the form ( |1 . 1 1|) , ( |1.13|) , 
which are finite for m = 0, see [31-34]. We shall always set m = and study the dependence 
of the scaling functions f(Mr) = f(Mr, mr = 0) on the only remaining argument Mr in the 
asymptotic region Mr << 1. This can be performed within the framework of the general 
solution fl2.21|) of the RG equations with the aid of additional methods, see Sec.[IV[ 



III. RENORMALIZATION AND CRITICAL DIMENSIONS OF COMPOSITE 

OPERATORS 



The quantities entering into the left hand sides of Eqs.( |l.ll| ), ( |1.13 ) are two-point cor 



relation functions of composite fields (composite operators in quantum-field terminology) 
rather than multipoint correlators of the primary fields. In what follows, we use the term 
"composite operator" for any local (unless stated to be otherwise) monomial or polynomial 
constructed from primary fields and their derivatives at a single point x = (t, x). Examples 
are 6 n (x), [d i d{x)d i d{x)] n , d i 6(x)d j 6(x), 9'(x)V t 9(x) and so on. 

Since the arguments of the fields coincide, correlation functions with these operators 
contain additional UV divergences, which are removed by additional renormalization pro- 
cedure, see e.g. IJ. For the renormalized correlation functions standard RG equations are 
obtained, which describe IR scaling with definite critical dimensions = A[F] of certain 
"basis" operators F. Owing to the renormalization, A[F] does not coincide in general with 
the naive sum of critical dimensions of the fields and derivatives entering into F. 

Detailed exposition of the renormalization of composite operators for the stochastic 



Navier-Stokes equation is given in the review paper [12|, see also [13-17,19-21], below we 
confine ourselves to only the necessary information. 

In general, composite operators are mixed in renormalization, i.e., an UV finite renor- 
malized operator F R has the form F R = F+ counterterms, where the contribution of the 
counterterms is a linear combination of F itself and, possibly, other unrenormalized opera- 
tors which "admix" to F in renormalization. 
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Let F = {F a } be a closed set, all of whose monomials mix only with each other in 
renormalization. The renormalization matrix Zp = {Z a p} and the matrix of anomalous 
dimensions 7^ = {70/3} for this set are given by 

F a = Y. Z <*pF* If = Z/V^Zp, (3.1) 

P 

and the corresponding matrix of critical dimensions A^ = {A a ^} is given by Eq. (|2.19j ), in 
which d F , d F , and dp are understood as the diagonal matrices of canonical dimensions of the 
operators in question (with the diagonal elements equal to sums of corresponding dimensions 
of all fields and derivatives constituting F) and 7* = 7(5"*) is the matrix (|3.1|) at the fixed 
point. 

Critical dimensions of the set F = {F a } are given by the eigenvalues of the matrix Ap. 
The "basis" operators that possess definite critical dimensions have the form 

= E UaflFf , (3-2) 
P 

where the matrix Up = {U a p} is such that A^ = UpApUp 1 is diagonal. 

In general, counterterms to a given operator F are determined by all possible 1- 
irreducible Green functions with one operator F and arbitrary number of primary fields, 
T = (F{x)<&{x\) . . . $(£2)). The total canonical dimension (formal index of divergence) for 
such functions is given by 

dr = dp — N$d<s>, (3.3) 

with the summation over all types of fields entering into the function. For superficially 
divergent diagrams, dr is a nonnegative integer, cf. Sec.[TT]. 

Let us consider operators of the form 6 n (x) with the canonical dimension dp = —n, 
entering into the structure functions ( 1.11 ). From Table | in Sec.|l| and Eq.( [3.3|) we obtain 



dr = —n + Ng — N v — (d + 1)^/, and from the analysis of the diagrams it follows that the 
total number of the fields 9 entering into the function V can never exceed the number of the 
fields 9 in the operator 9 n itself, i.e., Nq < n. Therefore, the divergence can only exist in the 
functions with N v = 0, Ngi = 0, and arbitrary value of n = Ng, for which the formal index 
vanishes, dr = 0. However, at least one of Ng external "tails" of the field 9 is attached to 
a vertex 9'{wd)9 (it is impossible to construct nontrivial, superficially divergent diagram of 
the desired type with all the external tails attached to the vertex F), at least one derivative 
d appears as an extra factor in the diagram, and, consequently, the real index of divergence 
is necessarily negative, see Sec.|J. 

This means that the operator 9 n requires no counterterms at all, i.e., it is in fact UV 
finite, 9 n = Z [9 n ] R with Z = 1. It then follows that the critical dimension of 9 n (x) is simply 
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given by the expression ( |2.19| ) with no correction from 7^ and is therefore reduced to the 
sum of the critical dimensions of the factors: 



A[6 n ] = nA[e] = n(-l + e/2). (3.4) 

We note that this relation was not clear a priori and it is a specific feature of the model 
Q2.1|). For example, in the standard model 4 of the theory of critical phenomena, the 



critical dimensions of the field <p and the composite operators <p and <p are completely 



independent, and they determine independent critical exponents i], u, and u, see e.g. |L0 
The relation analogous to (|3.4|) is valid for the powers of the velocity field of the stochastic 
Navier-Stokes equation, where it is a consequence of the Galilean symmetry of the model, 
see e.g. |TJ,|TJ,|lJ] . 



Now let us turn to the operators F n = [di6di6} n with oIf = 0, d F = —n. They enter into 
the left-hand sides of Eqs.( |l.l3| ) and, as we shall see in Sec.[TV], it is their critical dimensions 



that determine the anomalous exponents in ( |1 . 1 1| ) and ( |1.13| ). 

In this case, from Table | in Sec.[TI| and Eq.(|3.3|) we have dr = Ng — N v — (d + l)Ngi, 
with the necessary condition Ng < 2n, which follows from the structure of the diagrams. It 
is also obvious from the analysis of the diagrams that the counterterms to these operators 
can involve the fields 8, 6' only in the form of derivatives, d8, d6', and so the real index 
of divergence has the form d' r = d-p — Nq — Nqi = —N v — {d + 2)Ng>. It then follows 
that superficial divergences exist only in the Green functions with N v = Ngi = and any 
Ng < 2n, and the corresponding operator counterterms are reduced to the form Fk with 
k < n. Therefore, the operators F n can mix only with each other in renormalization, the 
corresponding infinite renormalization matrix Zp = {Z n k} is in fact triangular, Z n k = for 
k > n, and the critical dimensions associated with the operators F n are determined by the 
diagonal elements Z n = Z nn (in contradistinction with the case of operators 6 n , they are not 
equal to unity here). 

Explicit calculation of the constants Z n in the MS scheme in the two-loop approximation 
gives: 

n(d-l)(d + 2n) 
2e' d{d + 2) 
u 2 n{n-l)(d-lf(d + 2n)(d + 2n + 2) 
+ 8? d 2 {d + 2) 2 + 



u 2 n{n — l)(d — 1) 
+ 2e ' d 2 {d + 2) 2 {d + A) 



-(rf + 4)(rf + l)^ 



+3(rf - l)(d + 2n)h 1 (d)/A - id+1) ^^ n 2) h 2 (d) 



(d + 2) 

0(u 3 ), (3.5) 



where we have written u = gCd with the coefficient Cd from ( p.lOj ), and 
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/ii(d) = F(l,l,d/2 + 2,l/4), 
h 2 (d) = F(l, 1, d/2 + 3, 1/4) = /^(d + 2), 

where F(. . .) is the hypergeometric series (see e.g. [[46]]): 

„, , , , , ab a(a+l)b(b+l) z 2 

F(a,b,c,z) = l + - z +^ j-L^- J - ■- + .... 

c c(c+l) 2! 



(3.6) 



(3.7) 



From (13.61) and 



it follows: 



fc!r(d/2 + 2) 
^ 4 fc r(d/2 + 2 + A;) 



/12(d) = /ii(d + 2). 



(3.8) 



Representations ( [3.8D are convenient for numerical computations of the functions h\^(d). 
From ( |3.5| ) for the anomalous dimension r y n = T>^ In Z n in the two-loop order we obtain: 

unid — l)(d + 2n) u 2 n(n — 1) 

7n(#) = „ „ , . ^ + — X 



2d(d + 2) 



(d + 4)(1 - d 2 ) 3(d - l) 2 (d + 2n) , , (d 2 - l)(d + 3n - 2) 

H rm~ hi\d) — — — — — — ri2\d) 



d(d+2f 



4d(d + 2) 



d(d + 2)(d + 4) 



d(d+ 2)(d + 4) 

+ 0(w 3 ), (3.9) 



and for the corresponding critical dimension A„ at the fixed point ( 2.1 2| ) using ( [2.19 ) and 
d F [F n ] = 0, dp[F n ] = —n (see Table | in Sec.|l|) we have: 



A n = ne + 7* 



2n(n — l)e 



+ 



e 2 n(n — 1) 



d + 2 (d- l)(d + 2) 3 (d + 4) 2 



x 



-4(d + l)(d + 4) 2 + 3(d - l)(d + 2)(d + 4)(d + 2ri)h 1 (d) 



-4(d + l)(d + 2)(d + 3n - 2)h 2 {d) 



+ 0(e 3 ). 



(3.10) 



Expression ( |3.10| ) is simplified for any integer value of d owing to the fact that the series 
in (|3.7|) reduce to finite sums, see [fl(| : 



ht(d) = 2(d + 2) 
for any even value of d and 
/11(d) = 2(d + 2) 



k-2 



d/2+1 

-3)^ln(4/3) + g- /2 _ fc + 2 



-3) 



(3.11) 



(d+l)/2 /_o\(d+l)/2-A 
(_ 1 )(«*+l)/2. 3 d/2-l. 7r + 2 £ I <*J 

fe=l 



2k -1 



(3.12) 



for any odd value of d, which for d = 2 and d = 3 gives the results announced in ( |1.14|) , 

(EH)- 

We explain in Sec.[IV| that the critical dimensions A n from ( |3.10| ) are nothing else than 
the anomalous exponents entering into relations fll-11 ), (|1.13| ), and here we only note that 
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the first order in e of the expression (|3.10| ) coincides (up to the notation) with the result 
obtained in for n = 2 and in |34| for arbitrary n, and that the first term of the expansion 



in 1/d of the expression ( |3.10| ) coincides with the result ( |1.12| ) obtained in [pl] , p2| . We also 



note that the e 2 term of the expression (|3.10| ) behaves as 1/d 2 for 1/d — > and therefore 



gives no contribution to the first order of the 1/d expansion [this follows from the relation 
h,2(d) = h\(d + 2) = 1 + 0(l/d), which is obvious from ( |3.Sp l. This fact suggests that the e 
expansion for the dimension A n in the model fl2.1p is "better" than the 1/d expansion in the 
sense that a given order of the 1/d expansion is contained completely in the corresponding 
order of the e expansion, but not the reverse. We also note that the n dependence of the 
quantity A n in the second order of the e expansion is no longer reduced to the simple factor 
n(n — 1). 

The result Ai = in ( |3.10[ ) is in fact exact, in agreement with the exact solution for 



the two-point structure function obtained in [27]. Within the RG approach this can be 



demonstrated using the Schwinger equation of the form 

J V<$>5 \9{x) exp S R ($) + A$\ /56'(x) = (3.13) 

(in the general sense of the term, Schwinger equations are any relations stating that any 
functional integral of a total variational derivative is equal to zero, see e.g. Sec. 7 of [[47j[). 
In (|3.13|) , Sr is the renormalized analog of the action ( |2.1|) , and the notation introduced in 
( j2.2| ) is used. Eq. fl3.13|) can be rewritten in the form 

((9'D e 9 - V t [9 2 /2] + uZ u A[9 2 /2] - vZ v F x )) A = -A e ,8W R {A)/5A e . (3.14) 

Here Dg is the correlator (|1.8|) , ((...)) a denotes the averaging with the weight exp[S , ^($) + 
A&], Wr is determined by the Eq. ( |2.2| ) with the replacement S — > Sr, and the argument x 
common to all the quantities in ( |3.14|) is omitted. 

The quantity ((F)) a is the generating functional of the correlation functions with one 
operator F and any number of fields therefore the UV finiteness of the operator F 
is equivalent to the finiteness of the functional ((F)) a- The quantity in the right hand 
side of Eq. (|3.14j) is finite (a derivative of the renormalized functional with respect to finite 
argument), and so is the operator in the left hand side. Our operator Fx does not admix in 
renormalization to the operator 6'DgO (Fi contains too many fields 6), and to the operators 
Vt[0 2 /2] and A[6' 2 /2] (they have the form of total derivatives, and F± is not reduced to this 
form). On the other hand, the operator 6'Dg6 does not admix to Fx (it is nonlocal, and F\ is 
local), while the derivatives Vt[# 2 /2] and A[# 2 /2] do not admix to F\ owing to the fact that 
each field 9 enters in the counterterms of the operators F n only in the form of derivative 89 
(see above). Therefore, all three types of operators entering into the left hand side of the 
Eq.( |3.14D are independent, and they must be UV finite separately. 
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Since the operator vZ v F\ is UV finite, it coincides with its finite part, i.e., vZ v F\ = zaF/*, 
which along with the relation F\ = Z\F^ gives Z\ = Z^ 1 and therefore 71 = —7^. For the 
critical exponent A x = e + 7^ we then obtain Ai = exactly (we recall that 7* = e, see 
Sec.0). 

In the SNS model, critical dimensions of certain composite operators can sometimes 
be obtained exactly using various Schwinger equations and Ward identities for Galilean 
transformations, see [12-17] and [19-21]. In particular, exact critical dimension of the energy 
dissipation rate E(x) was found in [PO], see also [P1J] . The simple relation A[F ra ] = nA[E] 
for the powers of the dissipation rate was proposed in fUj . It was explained later in ||17|| , 
that this relation can not be considered reliable, and here we only note that even for the 
simple Obukhov-Kraichnan model, the critical dimension of E n is not a linear function in 
n, see ( |3.10| ), in contradistinction with the powers of the field itself, see ( |3.4j ) and p~^ , [13|JT5[ 
for the case of the SNS model. 

Critical dimensions of various tensor operators of the form di 1 9{x) . . . d,i n 6(x) can also be 
calculated in the second order of the e expansion from the same two-loop diagrams, which 



determine the constants (|3.5| ); only the symmetry coefficients differ from those for scalar 
operators. We shall confine ourselves to the second-rank irreducible traceless tensors of the 
form 



FZ^diPdjOFn-t-SijFn/d, 



(3.15) 



where F n are the scalar operators discussed above. The operators ( |3.15|) mix only with 
each other in renormalization, the corresponding renormalization matrix is triangular and 
its diagonal elements determine the corresponding critical dimensions A^. In the second 
order of the e expansion we have obtained: 



a: 



e[d{d + 1) - 2n(n - l)(d - 1)} 



e 2 n(n — 1) 



(d-l)(d + 2) {d-l)(d + 2f{d + Ay 

-4(d + l)(d + 4) 2 P + 3(d - 2){d + 4)(d + 2n)h 1 {d)P 1 - 



where we have written: 



-4(d + l)(d + 2){d + 3n - 2)h 2 (d)P 2 



0(e £ 



(3.16) 



Po 
Pi 
■ 1 - 



2n(n- l)(d- 1) 

d{d + l) 
n(d - l)(d + 2n)' 
d(2d- 1) 



2n(d- l)(d + 3n-2) ; 
and the functions hi^{d) are defined in ( |3.6|) . 



(3.17) 
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We note that the expression ( |3.16| ) for n — 1 coincides with the first two terms of the 



expansion in e of the exact result obtained in |H] (we explain below that the quantity 



A[ corresponds to 5 — 7 in the notation of pi|). For the case of the SNS model, critical 



dimensions of irreducible tensor operators are studied, for example, in [[16 . 

Since the critical dimensions of the operators 9 n and (d6) n have been found, we can use 
the general expression ( |2.21|) in the case of the structure functions ( |1.11|) and pair correlators 



The structure function S , 2n (r) = (\0(x) — 9(x')] 2n ) is represented as a sum of pair correla- 
tors (8(x) k 8(x') m ) with fixed value of k + m = 2n and equal canonical and critical dimensions 
d% = 2nd% = -n, A F = 2nA e = n(-2 + e). Then from Eq. (|2T21|) with m = it follows: 

S 2n {r)~D- n r n ^f n (Mr), (3.18) 

with certain, as yet unknown, scaling functions f n (Mr). 

The nth power of the dissipation rate is represented as a finite linear combination of 
basis operators (|3.2j ) with definite critical dimensions A& = ke + ^l given in ( p.lOj ), with the 
necessary condition k < n. Therefore, for the pair correlators ( |L13j ) of the powers of the 
dissipation rate we obtain from ( |2.21| ): 

n m 

(E n (x)E m (x>)) = ^ +m (F n (x)F m (x')) * E J2( Ar r Ak ~^h,i (Mr), (3.19) 

k=0 1=0 

with certain scaling functions fk,i (Mr) and the UV scale A defined in ( p.,10| ). The leading 
term of the asymptotic behaviour of the expression ( 3.19Q in the IR region Ar >> 1 is given 
by the contribution with minimal A*. + A;, i.e., 

(E n (x)E m (x>)) ~ (Ar)- A "- A ™/ n , m (Mr). (3.20) 



We note that the structure functions ( p.!8| ) and the correlator ( p. 20 ) for n = m = 1 are 



independent of the diffusivity coefficient, or, equivalently, of the UV scale A (we recall that 
Ai =0). This statement is an analog of the Second Kolmogorov hypothesis (independence 
of the viscosity coefficient in the inertial and energy-containing ranges) for the real fully 
developed turbulence, see e.g. M and [0]. Within the RG approach to the SNS model, the 



Second Kolmogorov hypothesis was established in fl3|l , see also |12|,|1 



Expressions like ( |3.18j ), ( |3.20| ) can easily be written down for any single-time pair cor- 



relator, provided its canonical and critical dimensions are known. Let us give two more 
examples. For the second-rank irreducible tensors E' n = vqF^ with the operators F' n given 
in ( |3.15|) we obtain (dropping the vector indices): 



(K{*)K>W)) = ^ +m (F' n (x)F' m (x')) * (Ar)- A "- A -/:, m (Mr), (3.21) 
and for the mixed correlator of the scalar E n and the tensor E' n we have: 
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(E n (x)E' m (x')) = u^{F n {x)F' m {x')) ~ (Ar)- A "- A ™/; m (Mr), (3.22) 

with the dimensions A n , A' n given in (|3.10| ), (|3.16| ) and certain scaling functions f nm (Mr), 
f" m (Mr) (the prime is not a derivative here). 



IV. OPERATOR PRODUCT EXPANSION AND ANOMALOUS SCALING 



From the viewpoint of the renormalization group, the expressions Q2.21 ) and ( 3.18|) — (|3.22|) 



for any functions /(Mr) correspond to IR scaling in the region L = M~ l , r >> I = A -1 
for arbitrary fixed value of Mr, with definite critical dimensions Ap. The inertial range 
I « r « L corresponds to the additional condition Mr << 1, and representations like 
( |1.11| ) and ( |1.13| ) should be understood as certain additional statements about the explicit 
form of the leading terms of the asymptotic behaviour for Mr — > 0. 

In the theory of critical phenomena, the asymptotic form of scaling functions for M — > 
is studied using the well known Wilson operator product expansion (OPE), see e.g. P,[10 



the analog of L = M 1 is there the correlation length r c . This technique is also applied to 



the theory of turbulence, see e.g. [12,13 



According to the OPE, the single-time product Fi(xi)F2(x2) of two renormalized oper- 
ators at x = (xi + x 2 )/2 = const , and r = x x — x 2 — > has the representation 

F 1 (x 1 )F 2 (x 2 ) = £C a (r)F*(x,t), (4.1) 

a 

in which the functions C a are the Wilson coefficients regular in M 2 and F^ are all possible 
renormalized local composite operators of the type Q3.2|) allowed by symmetry, with definite 
critical dimensions A a . 



The renormalized correlator (^(xi)^^)) is obtained by averaging (4.1) with the 
weight exp S R , the quantities (F£) = M da v d ™a a (g, M/fi) involving dimensionless functions 
a a (g,M/fi) appear on the right hand side. Their asymptotic behaviour for M/fi — > is 
found from the corresponding RG equations and has the form 

(^>ocM A ". (4.2) 

From the operator product expansion (|4.1|) we therefore find the following expression for the 
scaling function /(Mr) in the representation (|2.21|) for the correlator (-Fi(xi)F 2 (x2)): 



f{u) = Y,Mu)u A % u = Mr, (4.3) 

a 

with coefficients A a , which are regular in u 2 , generated by the Wilson coefficients C a in 
(|4.1|), which are regular in M 2 . 
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The leading contributions for u — > are those with the smallest dimension A a and in 
the e expansions they are those with the smallest d a = d[F a ] for e = 0. We shall term the 
operators with A a < 0, if they exist, dangerous [p~2| , p~3[1 , as they correspond to contributions 
to ( |4.3|) which diverge for u — > 0. 

In the standard model 4 of the theory of critical behaviour [TO], there are no problem 
of dangerous operators within the e expansions, because in that model A a = n a + 0(e), 
where n a > is the total number of fields and derivatives in F^. The operator F = 1 has 
the smallest value n a = 0, but it gives a contribution to ( |4.3| ) which is regular in u 2 and 
has a finite limit as u — > 0. The first nontrivial contribution is generated by the operator 
2 with n a = 2, it has the form u 2+ °^ and only determines correction to the leading term 
generated by the operator F = 1, which vanishes at u — > 0. 

We note that for a Galilean invariant product Fi(xi)Fz(xz), the right hand side of 
Eq. (|4. 1|) can involve any Galilean invariant operator, including tensor operators, whose in- 
dices would be contracted with the analogous indices of the coefficients C a . Without loss of 



generality, it can be assumed that the expansion is made in irreducible tensors (see Sec.|Tl 
for examples), so that only scalars contribute to the correlator (FiF 2 ) because the averages 
(Fa) for nonscalar irreducible tensors are zero. For the same reason, the contributions to 
the correlator from all operators of the form dF with external derivatives vanish owing to 
translational invariance. 

In our case, contributions from the operators like 9 n with negative dp are also forbidden 
by the invariance of the correlators ( |3.1S| )- p.22] ) with respect to the shift 9(x) — > 6 l (a;)+const 
of the field 9. 

The leading terms of the asymptotic behaviour of the scaling functions in ( 3.18Q — (|3.22|) 
for Mr — > are therefore determined by the scalar operators F n = [di9di9] n with the minimal 
canonical dimension dp = 0, see Sec. [II]. From the analysis of the diagrams it follows that the 
number of the fields 9 in the operator F n entering into the right hand sides of the expansions 
Q4.1D can never exceed the total number of the fields 9 in their left hand sides. Therefore, 
only finite number of operators F n contributes to each operator product expansion in the 
model (|2.1J ), and the asymptotic form of the scaling functions in (|3.18f )- (|3.22|) for u — > is 
given by the expression: 



A' 



f(u) = Y,C n u A " + ..., (4.4) 



n=0 



where iV is the total number of the fields 9 in the left hand sides, C n = C n (e, d) are numerical 
coefficients, and the dots stand for corrections of order u 2+0( - £ ^ vanishing as u — > 0. The 
leading term for u — > is determined by the operator with minimal A n , i.e., with maximal 
n (owing to the fact that the dimension A ra within the e expansion decreases monotonically 
with n, see ( |3.10 )). As a result, for the scaling function f n in the representation ( p . 1 8[ ) for 



21 



the structure function Son we obtain: 



f n ~ const u An for u —> 0, 



(4.5) 



and for the correlators (|3.20|)-(|3.22|) we have: 



fn,mj fn,mi fn 



const U 



for u — > 0, 



(4.6) 



with the critical dimensions A n given in ( 13.101 ). 

Therefore, we have derived expressions ( 1.11|) , (|1.13| ), and relate the corresponding 
anomalous exponents to the critical dimensions of the composite operators E n (x). 



V. DISCUSSION AND CONCLUSION 

We have shown that the renormalization group combined with the operator product 
expansion establishes the existence of anomalous scaling in the model ( |1.7| )-( |I~9|) for the 
advection of a passive scalar by a Gaussian velocity field and allows the corresponding 
anomalous exponents to be calculated in the form of series in e. 

The distinguishing feature of the model (1.7)- (|1.9|) that explains the origin of the anoma- 



lous scaling is the existence of dangerous composite operators with negative critical dimen- 
sions. They dominate the asymptotic behaviour of the scaling functions and lead to singular 
dependence of the correlation functions on the IR scale M for M — ► 0, in contrast to the 
standard models of critical phenomena, in which all the nontrivial operators have positive 
critical dimensions and only determine vanishing corrections to the leading finite contribu- 
tion from the simplest operator F — 1, see e. g. In contradistinction with the SNS 
model, the dangerous operators in the model (|i~^ )-(|l~^) occur already for asymptotically 
small values of e and only finite number of these operators contribute to operator product 
expansion for any given correlation function. 

The set of expressions fl3.18| )-( |3.22| ), (|4.4|) - ([4TB]) gives the complete description of the 



IR asymptotic behaviour of the Green functions in the model (p..7|) - ([L.9|) : Eqs. (|3.18|) - (|3.22|) 



describe the asymptotic form of the structure functions and pair correlators in the IR region 
L,r » I = A and determine their dependence on the UV scale A, while Eqs.( |4.4|) -( [4.6| ) 
give the asymptotic form of the corresponding scaling functions upon additional restriction 
Mr << 1 and determine the dependence on the IR scale L = M~ l . All the critical di- 
mensions A n , entering into (|3.18|) - (|3.22|) , ([4.4|) - ([4.6|) have been calculated in the second 



order of the e expansion, see fl3.10p , ( |3.16| ), while the critical dimensions of the model pa- 
rameters, primary fields, and their powers have been found exactly, see ( [2.18|) , (|2.19|) , ( [2.20|) 
and Q. 
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It should be stressed that the asymptotic expressions Q4.5|) , ( |4.6|) result from the fact 
that the critical dimensions A n are negative and that |A n | increases monotonically with 
n. This is obviously so within the e expansion, in which the sign and the n dependence 
of the dimensions are determined by the first-order terms, while the higher-order terms are 
treated as small corrections. However, for finite values of e the higher-order terms can, 
in principle, change these features of the dimensions. Indeed, the n 3 contribution in the 
second-order approximation for A n is positive, see e.g. ( |1.14j ), ( |1.15| ), and so A n becomes 
positive, provided n is large enough. Of course, this conclusion is based on the second-order 
approximation of the e expansion and is therefore not reliable: higher-order terms of the e 
expansion contain additional powers of n, and the correct analysis of the large n behaviour 
of the dimensions A n requires resummation of the e series with the additional condition that 
ne ~ 1, but we know of no model in which such a resummation has been performed. 

The comparison of the Eqs. ([3.18|) -( [3.22|) , ( |4.5| ), ( 4~6|) with the corresponding expressions 
in [31-34] shows that A n = — p 2n in the notation of [3~3|, A 2 = —A, A[ = 7 — S in the 



notation of [31|, and A n = — A2 n in the notation of |32]; our results for A„, A^ are in 
agreement with the results obtained in [31-34] for the structure functions S < 2 n (r) and the 
pair correlators ( |3.20| ), and in |JT] for the correlators ([3.2 1|) , ( |3.22j ) with n = m = 1 within 
the first order of the expansions in e and 1/d. 

It is noteworthy that the set of scalar operators F n = [di6di6} n is "closed with respect to 
the fusion" in the sense that the leading term in the OPE for the pair correlator (F n F m ) is 
given by the operator F n+m from the same family with the summed index n + m, see Q3.20D , 
(fl.6|). This fact along with the inequality A„ + A m > A„ +m , which is obvious from the 
explicit expressions for A n , can be interpreted as the statement that the correlations of the 
local dissipation rate in the model (|1.7| )- (|1.9|) exhibit multifractal behaviour, see j|9| and 
We note that the same relation ensures the fulfillment of the Holder inequality for the 



structure functions ( |1.11| ). 



An important question is that of the universality of anomalous exponents, see e.g. 35 



It is clear from the RG analysis, that the exponents A n do not depend on the choice of the 
correlator ( |1.9|) (this correlator does not enter at all into the UV divergent diagrams which 
determine renormalization constants), and that they are insensitive to the specific form of 
the IR regularization in the correlator (|1.8| ) (renormalization constants do not depend on the 
choice of the IR regularization, see e.g. ||]). However, the anomalous exponents can change 
if the function S(t — t') in the correlator of the velocity field is replaced by some function 
with finite width, i.e., the velocity field has small but finite correlation time [|36|. The RG 
approach to this problem will be presented elsewhere, and here we only mention another 
possible modification of the model ( |1.7|) -( fL9|) , that of a "slow" velocity field. In this case, 
the function 5(t — t') in ( |1.9| ) is replaced by the unity, so that the velocity correlator is time 
independent. 



23 



The RG analysis can be directly extended to this model to prove that its Green functions 
also exhibit anomalous scaling behaviour and the corresponding anomalous exponents can 
be calculated in the form of series in e = e + 2. The critical dimensions A^ = 2 — e/2, 
Aq = — 1 + e/4 are found exactly, and for the structure functions defined in ( |1 . 1 1| ) we have 
obtained: 

5 2n ~Do" /2 r n(2 " l/2) (Mf) A ", (5.1) 

where 

A n = -en(n- l)/(d + 2) + 0(e~ 2 ). (5.2) 

The expressions (|1.13|) remain valid with the replacement A„ — > A n . We note that the 
velocity field with the dimension A[v] = 1 — e/2 becomes dangerous for e > 2 (which 
corresponds to the IR divergence of the integral in Q1.SQ with m — 0), and so become 
all its powers. However, these operators are not Galilean invariant and therefore give no 
contribution to the operator product expansions of the structure functions and correlators 
d). 
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TABLES 



TABLE I. Canonical dimensions of the fields and parameters in the model (2.1). 
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